"'--■ft 


HOLOGRAPHIC  INTERFEROMETRIC  SURVEY  OF  BOUNDARY 
LAYER  TRANSITION  IN  AN  AXISYMMETRIC  FREE 
CONVECTION  THERMAL  BOUNDARY  LAYER 


THESIS 


AFIT/GAE/AA/79D-12 


Gregory  S.  Meserve 
2nd  Lt  USAF 


D  C 

nn 

,  .  - i- - :.■> - l 

I  %  JAN  29  1980 


(  i 


Li  Li  L 


:-U  TTE 

A 


Approved  for  public  release;  distribution  unlimited. 


k 


ORIGINAL  CONT/L  -  PLATES:  ALL  DO<r 

REPRODUCTION  >  i..L  Lt.  iN  BLACK  AND  WHfTf 

. — - - - - — — - - — - - - - — — - 


AFIT/GAE/AA/79D-12 


Preface 


At  first  the  main  objective  of  this  study  was  to  train  its  author  in 
the  field  of  pulse  laser  holography.  Upon  graduation  the  Air  Force  Flight 
Dynamics  Lab  would  take  on  the  student  to  assist  on  projects  that  required 
optical  measurements. 

A.  G.  Havener  presented  the  objective  to  J.  Hitchcock,  professor  at 
the  AFIT  School  of  Engineering.  A  meeting  of  a  student  and  the  two  men 
was  established  and  initiated  the  beginnings  of  this  study. 

A  heated  vertical  cylinder  with  free  convection  conditions  was  chosen 
to  be  the  test  subject  because  experimental  data  is  not  known  to  exist. 

The  results  would  confirm  the  numerical  theoretical  solutions  that  were 
available. 

The  selection  of  the  vertical  cylinder  provided  interesting  and 
informative  conclusions,  but  the  main  objective  was  never  achieved  however. 
The  task  of  mastering  the  pulse  laser  holography  equipment  was  successful 
but  the  author  will  be  stationed  elsewhere. 

I  am  indebted  to  A.  G.  Havener  and  J.  Hitchcock  for  their  time, 
patience  and  insight.  The  experience  has  enhanced  my  awareness  and 
knowledge . 

Finally,  I  wish  to  thank  my  wife  Karen  who  has  witnessed  many  evenings 
alone.  Our  love  strengthens  as  our  list  of  life's  experiences  continues 
to  grow. 
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Abstract 

Unique  optical  data  for  the  laminar  and  transition  regions  of  an 
axisymmetric  free  convection  thermal  boundary  layer  are  presented.  The 
boundary  layer  encompasses  a  vertically  oriented  right  circular  cylinder 
kept  at  a  constant  wall  temperature  and  suspended  in  a  quiescent  air 
environment. 

Pulse  laser  holographic  interferometry  was  used  to  obtain  a  measure¬ 
ment  of  the  temperature  distribution  through  the  boundary  layer.  Schlieren 
methods  were  used  to  photograph  the  boundary  layer  to  provide  a  qualitative 
description  of  the  transition  region.  The  eddies  of  transition  flow  were 
the  important  aspect  observed  in  the  schlieren  measurements. 

An  approach  is  outlined  to  find  an  approximate  numerical  solution  for 
the  temperature  distribution  through  the  laminar  boundary  layer.  Integral 
equations  are  developed  for  a  differential  control  volume  in  the  boundary 
layer  after  which  assumed  forms  of  the  velocity  and  temperature  distribu¬ 
tions  are  used  to  reduce  the  integral  equations  to  first  order  non-linear 
differential  equations.  The  differential  equations  are  non-dimensionalized 
and  numerical  techniques  were  used  in  an  attempt  to  obtain  a  solution. 

Quantitative  data  from  the  laminar  region  resulted  in  a  discrepancy 
of  34.7  percent  between  theoretical  and  experimental  evaluation  of  a  local 
Nusselt  number.  Qualitative  observation  of  the  transition  region  reveals 
a  fringe  reversal  phenomena  that  could  possibly  be  associated  with  the 

fundamental  mechanics  of  transition.  Furthermore/ vortices  or  eddies  were 

A 

observed  in  the  transition  region. 


HOLOGRAPHIC  INTERFEROMETRIC  SURVEY  OF  BOUNDARY 
LAYER  TRANSITION  IN  AN  AXISYMMETRIC  FREE 
CONVECTION  THERMAL  BOUNDARY  LAYER 


L 

I.  Introduction 

Background 

To  date,  the  analyses  of  boundary  layers  for  vertical  cylinders  are 
numerical  and  are  essentially  unsupported  by  experimental  data.  Two 
studies  that  pertain  directly  to  this  effort  are  good  examples.  In 
1956,  Sparrow  and  Gregg  (Ref  9)  developed  a  numerical  solution  for  the 
ratio  of  vertical  cylinder  to  flat  plate  Nusselt  numbers  as  a  function  of 
axial  position,  and  in  1974,  Cebeci  (Ref  4),  solved  the  governing  system 
of  equations  for  laminar  free  convection  from  the  outer  surface  of  a 
vertical  cylinder  by  using  a  two-point  finite-difference  numerical 
technqiue.  Unfortunately  no  experimental  data  are  available  for  com¬ 
parison  and  conformation  of  these  numerical  solutions. 

The  lack  of  experimental  data  of  sensitive  flows  is  due  primarily  to 
a  lack  of  adequate  instrumentation.  Until  the  development  of  coherent 
optical  measuring  techniques,  probing  devices  had  to  be  used  to  obtain 
quantitative  field  data.  When  the  flow  field  is  sensitive,  the  intrusion 
of  a  probe  always  makes  the  measurements  suspicious.  Optical  techniques 
overcome  this  difficulty  because  the  beam  of  laser  radiation  emitted  from 
a  relatively  low  energy  laser  (less  than  one  joule,  typically)  is  a  non- 
disturbing  probe.  The  techniques  for  producing  holograms  developed  shortly 
after  the  Invention  of  the  laser,  and  since  then  steady  progress  has  been 
made  by  researchers  using  holography  in  conjunction  with  interferometry 
(  to  explore  many  different  areas.  In  gas  dynamics,  for  example.  Havener 
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(Ref  10)  and  Trolinger  (Ref  17)  have  successfully  developed,  perfected  and 
applied  holography  to  'imerous  wind  tunnel  investigations  during  the  past 
twelve  years. 
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The  objective  of  this  study  was  to  measure  temperature  distributions 
within  the  free  convection  thermal  boundary  layer  encompassing  a  vertically 
oriented  circular  cylinder.  Pulse  laser  holography  was  used  to  produce 
holograms  from  which  temperatures  and  Nusselt  numbers  were  extracted.  This 
information  was  compared  to  Sparrow’s  and  Gregg’s  numerical  analysis  to 
see  how  the  experimental  data  agree  with  the  numerical  data. 

Following  an  examination  of  the  laminar  boundary  layer,  new  and  unique 
data  of  the  transition  region  were  obtained.  These  data  we',  c  used  in  hope 
of  getting  a  better  understanding  of  the  transition  process.  Also, 
schlieren  photography  was  used  to  examine  qualitatively  the  transition 
boundary  layer.  Finally,  following  an  integral  approach  suggested  by 
Eckert  (Ref  7),  a  numerical  scheme  was  devised  to  predict  the  temperature 
distributions,  but  unfortunately,  satisfactory  results  of  this  approach 
aTe  unavailable  at  this  time. 


II.  Equipment 


The  vertical  cylinder  was  a  polished  copper  tube  six  feet  in  length 
with  a  one  and  five  eights  inches  outer  diameter.  The  cylinder  wall  was 
heated  uniformly  by  five  internally  fixed  quartz  infrared  lamps  connected 
in  series.  Each  lamp  was  separated  by  ceramic  spacers  which  give  added 
support  to  the  heating  elements  and  also  minimize  non-uniform  convection 
heating  inside  the  cylinder.  A  220  volt  variac  was  used  to  control  the 
cylinder  wall  temperature  by  regulating  the  voltage  supplied  to  the  lamps. 
Skematic  diagrams  of  the  cylinder  and  a  ceramic  spacer  are  shown  in  Figs 
1  and  2,  respectively. 

Thermocouples  of  chromel  -  alumel  are  imbedded  in  copper  rods  inserted 
in  five  of  the  ceramic  spacers.  The  rods  were  spring  loaded  and  have 
rounded  polished  contact  points  that  help  provide  good  thermal  contact 
between  the  copper  rods  and  the  inside  cylinder  wall.  The  thermocouples 
are  positioned  circumferentially  at  45  degree  spacings  around  the  cylinder 
axis,  and  an  ice  bath  is  used  for  the  reference  junction. 

A  Hi-Cam  high  speed  motion  picture  camera  was  used  in  conjunction 
with  a  regular  Toepler  schlieren  system  to  obtain  a  history  of  the  flow 
in  the  transition  and  turbulent  regions.  The  camera  was  set  to  expose  five 
hundred  frames  per  second  and  a  tricolor  filter  was  used  in  place  of  the 
schlieren  knife  edge  normally  used  in  schlieren  systems.  The  filter  was 
composed  of  blue,  green  and  red  transparent  plastic  strips  (see  Fig  3) . 

The  width  of  the  green  strip  was  approximately  the  same  as  the  width  of 
the  image  of  the  point  light  source.  In  use,  a  light  ray  refracted  off 
line  will  pass  through  the  red  or  blue  strips  depending  on  the  direction 
of  refraction  and  the  respective  orientation  of  the  filter.  For  this  test, 
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TO  VARIABLE  VOLTAGE 


Fig.  1  Internal  View  of  Cylinder 


Pig.  3  Composition  of  Color  Filter 


the  tricolor  filter  was  oriented  vertically  so  that  rays  traversing  the 
thermal  boundary  layer  pass  through  either  red  and  blue  strips  while  rays 
traversing  the  ambient  air  pass  through  the  green  strip.  Then,  when  the 
cylinder  was  photographed  with  white  light,  the  ambient  region  appears 
green  while  the  boundary  layer  appears  red  on  one  side  and  blue  on  the 
other.  Figure  4  is  a  typical  example  of  this  result. 

The  holographic  system  used  in  this  study  was  powered  by  a  pulse  ruby 
laser.  This  system  has  been  in  use  for  more  than  ten  years  and  is  docu¬ 
mented  in  several  references  (Ref  10) . 

The  cylinder  was  suspended  in  the  test  cabin  of  a  hypersonic  wind 
tunnel.  The  test  cabin  can  be  sealed  tightly  thereby  minimizing 
environmental  disturbances  which  can  easily  disrupt  the  free  convection 


thermal  boundary  layer.  Further,  this  test  cabin  has  eighteen  inch 
diameter  windows  which  permit^ large  aperture  viewing;  for  this  study. 


it  was  possible  to  view  eighteen  inches  of  cylinder  at  one  time.  The 
complete  system  schematic  diagram  is  shown  in  Fig  5. 


Fig.  4  Print  of  Cylinder  and  Boundary  Layer  with  Colored 
Filter  In  Use 


Pig.  5  Arrangement  of  Optical  Components  Used  to 
Produce  Holograms 
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III.  Theory  for  Data  Collection 

Interferometry  is  an  optical  process  that  can  be  used  to  gather 
experimental  data  on  refractive  index  variations  in  a  test  speciman  and 
holography  is  the  process  that  can  be  used  to  record  the  three  fundamental 
properties  of  light:  frequency,  intensity  and  phase.  The  hologram  is  a 
high  resolution  photographic  emulsion  in  which  the  full  content  of  a  light 
wave  reflected  from  or  transmitted  past  a  subject  can  be  stored.  In 
practice,  the  following  general  process  is  used  to  produce  a  hologram. 

Light  from  a  laser  is  divided  into  two  separate  beams.  One  beam  illuminates 
the  subject  which,  in  turn,  reflects  or  transmits  the  light  to  a  special 
film.  During  the  same  time,  the  other  beam  is  used  as  an  optical  reference 
which  directly  illuminates  the  same  film.  After  the  film  is  exposed, 
chemically  processed  and  reilluminated  with  only  the  reference  beam,  a 
three  dimensional  image  of  the  subject  can  be  seen. 

In  dual  hologram  interferometry,  a  hologram  of  the  ambient  condition 
and  one  hologram  of  the  test  condition  are  simultaneously  illuminated 
with  the  reference  beam,  thereby  producing  simultaneous  reconstruction 
of  the  object  waves.  The  reconstructed  ambient  waves  serve  as  a  reference 
against  which  the  reconstructed  test  waves  are  compared  in  the  formation 
of  the  interferogram. 

A  shift  in  the  phase  of  a  light  wave  is  revealed  in  an  interferogram 
as  a  fringe  shift.  A  fringe  is  a  line  of  constant  optical  path  length 
and  a  displacement  of  a  fringe  or  fringe  shift  is  caused  by  changes  in 
the  relative  phase  of  the  light  wave.  The  phase  of  the  wave  varies  when 
the  wave  traverses  a  volume  with  a  variable  refractive  index  field.  For 
gases,  refractive  index  variations  are  proportional  to  density  variations 
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and  consequently.  Interferometry  can  be  used  to  measure  density 
variations  in  gas,es.  The  governing  equation  of  this  process  is 


S 


i 


f(x,y,z)  d l± 


(1) 


where 

f(x,y,z)  is  the  density  function  —  1 

Pref 


i  is  a  subscript  denoting  an  arbitrary  straight 

line  through  the  density  field 

is  the  optically  measured  fringe  shift 
corresponding  to  the  i-th  line 

■t  is  the  geometric  length  of  the  i-th  line  through 

the  density  field 

G  is  an  optical  constant  prescribed  by  the  wavelength 

of  the  light,  by  the  composition  of  the  density 
field,  and  by  the  reference  density 

<G--F  '■ref) 

Equation  (1)  is  used  to  calculate  local  density  in  the  boundary  layer  from 
interferometric  fringe  shifts  obtained  from  the  holograms. 

Analytic  solutions  of  Eqn  (1)  are  possible  only  when  the  density 
field  is  two-dimensional  or  axisymmetric  which  is  an  important  assumption 
of  this  investigation.  To  obtain  the  density,  Eqn  (1)  must  be  inverted. 
Analytical  inversion  is  possible  when  the  number  of  independent  variables 
In  the  integrand  are  reduced  to  one.  Dependence  on  one  of  the  variables 
is  suppressed  by  representing  the  density  field  as  the  summation  of  thin 
slices  (for  this  case  horizontal  slices  perpendicular  to  the  cylinder  axis) 
within  each  slice  the  density  depends  only  upon  two  variables.  A  complete 


II 


description  is  obtained  by  combining  the  results  of  all  the  slices. 

Using  rectangular  cartesian  coordinates  and  letting  the  x  axis 
coincide  with  the  axis  of  the  cylinder;  the  density  field  is  sliced  into 
thin  planes  in  the  x-direction.  In  each  plane  the  density  depends  only 
upon  y  and  z,  The  i-th  line  through  the  field  can  be  expressed  in  polar 
coordinates  where  r  is  the  distance  from  the  origin  (axis  of  cylinder 
for  this  study)  and  0  is  the  corresponding  view  angle  (Fig  6) .  Equation  (1) 
can  be  written  for  a  slice  as 


S’(r.e)  -  G 


f(y,z)  dl± 


For  an  axisymmetric  density  field  all  views  of  the  field  are  the  same; 
therefore,  dependence  on  6  is  eleiminated  and  Eqn  (2)  becomes 


S(r)  =  G  Jt  f (y,z)  d l± 


where 


/  2  .  2 
y  +  z 


To  obtain  an  analytical  solution  to  Eqn  (3)  the  dependence  on  either  y  or 
z  must  be  eliminated.  To  accomplish  this,  the  density  field  is  assumed 
to  be  refractionless,  which  means-  the  bending  of  a  light  ray  from  a 
straight  line  through  the  density  field  can  be  neglected  because  this 
deviation  is  assumed  to  b.e  small  compared  to  the  width  of  the  field. 

i 

Taking  R  as  the  distance  from  the  origin  to  the  edge  of  the  boundary  layer 
and  0  equal  to  zero  so  that  r  is  equal  to  y,  Eqn  (3)  becomes 


S(y)  -  G 


fCr)  dr 


Pig.  6  General  View  of  Density  Field  For  Constant 
x-Plane 


Here,  y  is  constant  meaning  that  the  integration  is  carried  out  in  the 
2  direction  only.  By  using  the  transformations 


and  Eqn  (8)  becomes 


Where  the  coefficients  can  be  determined  from  Eqn  (14)  by  using  a  least 
squares  curve  fitting  routine  to  represent  the  measured  S(£)  values. 

For  this  study  the  pressure  is  assumed  to  be  constant  across  the 
boundary  layer  and  the  air  is  assumed  to  behave  as  an  ideal  gas.  Hence, 
the  temperature  through  the  boundary  layer  is  inversely  proportional  to 
the  density.  When  the  distribution  of  density  through  the  boundary  layer 
is  known,  the  temperature  distribution  can  then  be  determined. 
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IV.  Development  of  Mathematical  Theory 


Equation  of  Motion 

In  Fig  7  the  forces  and  momentum  flux  terms  that  act  on  a  differential 
control  volume  are  shown.  Here 


dw  =  weight  of  fluid  in  c.v.  =  g/gc 


prdrdx 


(MV)±n  =  momentum  flux  in 


=  2t r  J 

J  r0 

£ 


pv  ^rdr 
x 


2  5  I  2 

(MV)  out=  momentum  flux  out  =  2tt |  pv^  rdr  +  2ir  I  pvx  rdrdx 


dFr  =  shear  force  *  2tt  r  dx 

w  o 


PAi  =  pressure  force  acting  on  surface  at  x  ==  2ttP j  rdr 

r0 


*  2ttp/* 
J  Tr 


-  pressure  force  acting  on  surface  at  x  +  Ax 


2irP 


/lIdr  +  2^[P/>] 


dx 


The  equation  of  motion  which  is  derived  from  Newton's  Second  Law 
states  that  the  rate  of  change  of  momentum  of  the  fluid  passing  through 
the  control  volume  boundary  is  equal  to  the  sum  of  the  forces  acting  on 
the  control  volume.  With  regard  to  normal  boundary  layer  assumptions 
(  ’>>  )  the  integral  equation  of  motion  is 

(18) 


x  r 

*  J  rt 


pv  r  dr 
K  x 


-wR 

J  T, 


r  dr  -  g/gc  /  pr  dr  -  tw  t£ 

o  J  To 
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For  Newtonial  fluids  the  shear  stress  at  the  wall  is 


V  ( 


3v£ 

3r 


(19) 


The  pressure  gradient  is  determined  by  considering  a  fluid  element 
outside  the  boundary  layer.  The  pressure  terms  effecting  this  differential 
fluid  element  are  shown  in  Fig  8.  Evaluation  of  the  pressure  terms  yields 

g 

^  (P2irr  dr)  dx  =  2irr  dr  dx  (20) 


and 


9P 

3x 


“  8/gc  P, 


(21) 


Substituting  Eqns  (19)  and  (21)  into  Eqn  (18)  yields 


(P»  "  P>  rdr  +  ^ro  <  f?*  > 


(22) 


From  thermodynamics  the  coefficient  of  volumetric  expansion  is 
defined  as 

B  -  -  <“> 


and  for  small  temperature  differences,  8p  -  constant.  If  the  fluid 
properties  are  evaluated  at  a  mean  boundary  layer  temperature 

*  *j(t w  +  ta)  the  assumption  $p  -  constant  is  better.  Integrating, 

Eqn  (23) 


P„  -  P  ■  Mt  -  tj 


(24) 


17 


Pig.  7  Arisymmetric  Differential  Control  Volume  For 

Momentum  Integral  Equation  Derivation,  Steady  Plow 


P2TTrdr  4  (P2ITrdr)  rix 


j: 


dx 


1 

I 


dr 


_J 


t  r  g/Bc  Poo21h  drdx 

P2Hrdr 


Pig.  8  Differential  Fluid  Element  Outside  Boundary 


Assuming  the  flow  Is  incompressible,  end  substituting  Eqn  (24)  into 


Eqn  (22)  leads  to 


tif  f 

J  r  J  r 


(t  -  tj 


xdr  -  v  (f^*) 


(25) 


w 


By  defining  y  =  r  -  tq  the  coordinate  system  is  transferred  to  the 
surface  of  the  cylinder,  and  by  integrating  from  the  wall  to  the  edge 
of  the  boundary  layer  Eqn  (25)  becomes 


1_ 

r 

o 


(y  +  *o)  dy  = 


(26) 


Note  that  if  y/rQ  <  .1  Eqn  (27)  reduces  to  the  vertical  flat  plate 
equation  by  letting  y  -  0. 
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Energy  Equation 

For  the  differential  control  volume  shown  in  Fig  9,  the  contributing 
energy  terms  and  energy  flux  terms  are: 


dQ^n  =  heat  transfer  rate  =  q^  2Trrodx 


I.  =  rate  at  which  enthalpy  is  convected  into  c.v. 
in 


=  2ir  /  pvxirdr 


I  ^  =  rate  at  which  enthalpy  is  convected  out  of  c.v. 
out 


=  2it /  v;.irdr  +  2n  -r^-  /  pv  irdrdx  +  p  v  i 

7  x  9x  J  K  x  eet 


2nRdx 


The  kinetic  and  potential  energy  terms  of  the  mass  flow  into  and 
out  of  the  control  volume  are  neglected  assuming  that  their  contribution 
is  small  compared  to  the  enthalpy.  Furthermore,  steady  state  is  assumed 
and  no  work  is  being  done. 

The  energy  equation  in  integral  form  for  this  control  volume  is 


R 

—  —  f 

r„  3*  J 


PVrdr  +  T  '  "w 

o 


where  the  heat  flux  at  the  wall  is  given  by  Fourier's  heat  conduction 


law  as 


,  ,  3t  . 

S,  "  -k  <  3?  >w 
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Fig.  9  Axis ymme trio  Differential  Control  Volume  For 
Energy  Transfers  In  A  Steady  State  Flow 


Front  the  law  of  conservation  of  mass  applicable  to  this  differential 


control  volume , 


p  v  2ir  R  dx  = 
e  e 


■tf. 


pv  2rrr  dr  dx 
x 


A±  [ 

R  8x  / 

J  r 


pv  r  dr 
x 


Substituting  Fourier’s  heat  conduction  law  and  Eqn  (29)  into 


Eqn  (27)  yields. 


—  —  f 

To  3*  / 

J  r. 


pv  (i-i  )  r  dr  «  *k  (  7-  ) 

X  6  OX 


For  an  ideal  gas,  enthalpy  is  a  function  of  temperature  only  and 


can  be  expressed  as 


di  «=  c  dt 
P 


Assuming  c^  is  approximately  constant  over  the  temperature  range  in 
the  boundary  layer,  integration  of  Equation  (31)  leads  to 


i-i  ”  c(t  -  t  ) 

00  p  00 


where  i  «  i  and  the  thermal  diffusivity  a  is,  a  =  k /pc  ,  Eqn  (30) 

»  e  p 


can  be  written  as 


ii  [ 
ro  J 


vx  (t  ~  **>r  dr  *  "  a 
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Using  the  transformation  y  =  r-ro  as  before,  integration  of  Eqn  (33) 
between  the  wall  and  the  edge  of  the  boundary  layer  leads  to 

t- 

r 5 

roh  I  &  +  v-iy  -  -  «0w  (34) 

J  o 

Similar  to  the  equation  of  motion,  note  that  if  y/rQ  1  .1,  the 
energy  equation  reduces  to  the  vertical  flat  plate  equation  by  letting 
y  -  0. 

Equations  (26)  and  (34)  are  a  pair  of  coupled  integral  equations 
derived  for  an  axisymmetric ,  differential  control  volume.  From  here,  the 
objective  is  to  transform  these  integral  equations  into  non-linear  first 
order  differential  equations.  Then,  the  non-linear  first  order  differential 
equations  are  non-dimensionalized.  Finally,  numerical  approximate 
techniques  are  used  at  an  attempt  to  solve  the  two  non-dimensional,  non¬ 
linear,  first  order  equations  simultaneously. 

In  a  theoretical  development  for  a  vertical  oriented  flat  plate 
free  convection  boundary  layer  analysis  which  is  similar  to  the  development 
here,  Eckert  assumed  the  form  of  the  velocity  and  temperature  profiles 
to  evaluate  the  integrals  in  the  pair  of  coupled  flat  plate  integral 
equations.  The  profiles  assumed  by  Eckert  will  be  used  here.  Accordingly, 
the  velocity  profile  is 

Vx  "  y/5)^  (35) 

with  v  being  an  arbitary  function  with  the  dimensions  of  velocity.  The 
maximum  value  of  velocity  v^  occurs  when  y  ”  5/3.  At  this  distance. 
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(v  )  *  4/27  v.  The  boundary  conditions  of  v  -  0  at  y  ■  0  and  y  *  6 

x  max  x  *  * 

are  also  satisfied.  The  temperature  profile  is  approximated  by  the 
parabola 

t  -  =  (tw  -  tj  a  -  y/6)2  (36) 


and  the  boundary  conditions  t  =*  t  at  y  =  0  and  t  ®  t^  at  y  =  6  are  also 

satisfied.  These  profiles  are  shown  in  Fig  10. 

Let  8  =  t  —  t  and  0  “  t  -  t  .  Substitution  of  the  velocity  and 

00  w  w  00 

temperature  profiles  into  Eqns  (26)  and  (.34)  will  result  in  the  follow¬ 
ing  first  order  differential  equations.  Note  that  the  partial  derivatives 
with  respect  to  x  can  be  written  as  total  derivatives  since  the  variables 
are  functions  of  x  alone.  The  resulting  equation  of  motion  is 


1  d_ 
280  dx 


2  .,5  ,  8.1  QQ  f5 

v  6(—  +  3)J«  g00w 


"Arf 


and  the  energy  equation  is 

v  - 


1  _d_ 

105  dx 


2a 

6 


(37) 


(38) 


Using  the  non-dimensionalizing  factors  prescribed  by  Sparrow  and 
Gregg, 


x 


+  VX 

v  "*5? 


+ 

X 


X 

r0Grx*s 


where  the  Grashof  number,  Gr,  is 


CYLINDER  WALL 


BOUNDARY  LAYER  THICKNESS 


.Fig.  10  Assumed  Temperature  and  Velocity  Profiles 


Equations  (37)  and  (38)  are,  respectively. 


1  d 

1+3 


x+5  <5+  v^d  +  |  6+  x+)  j 


140  6+x+(l+l/4  6+x+)-420 

6 


and 


(39) 


1  d 
***  d? 


[*+3  6+v+  a  + 1 SV>] 


(40) 


Upon  completion  of  the  differentiation,  Eqn  (39)  is 

xVv^U  +  f  z+5+)  +  2x+«+Vu  +  i  x+«+)  |4  - 


dx'* 


1403+2  (1  +  |  x+6+)  -  420  v+  -  S6+2  v+2  (1  +  x+d’1')  (41) 


and  Eqn  (40)  is 

«w  (1  +  4  xV)  +  xV2  (1  +  4  xV)  44 

7  ix+  7  <1*+ 


240 

Pr 


-  3«+2v+(1+~  x+6+) 


(42) 


By  defining  the  following  coefficients 


Cx  -  x+6  V2  (1+|  xV) 


C2  -  2 x+6+2v+  <1+|  x+d+) 

C3  -  1405+2  (1  +  i  x+d+)  -  420  v+-5d+2v+2(l+^j-x+d+) 


C4  -  x+6+v+  (1  +  ^  xV) 
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_  +.+^  /i  ,  2  +.+* 

C5«xoi  (1  +  y  x  6  ) 

_  240  t  8  +*+  ,  1  \ 

C6  "  TT  “  36  V  (  2l  x  6  +1) 


the  equation  of  motion  and  the  energy  equation  can  be  written,  respectively, 
as 


c  +  c  K 

1  dx+  2  dx+ 


(43) 


d6 

d? 


+  C 


dv 


5  dx^ 


=  C. 


(44) 


-f  +  + 

Initial  values  for  6  and  v  can  be  calculated  by  setting  x  =  0  in 

the  first  order  differential  equation  of  motion  and  equation  of  energy  . 

The  two  equations  can  be  treated  as  two  simultaneous  equations  with  the 

two  unknowns  being  6  and  v  .  The  resulting  equations  for  the  initial 

conditions  are 


5+  =  (240)**  (21  +  Pr)i* 

*  O 


(45) 


•  + 


(SO)* 


(|?  +  Pr)1* 


(46) 


The  Initial  slopes,  i.e.  (d<5  /dx  )^  and  (dv  /dx  )^  can  be  calculated 
by  evaluating  the  limit  as  x+  approaches  zero  in  the  two  non-dimensional 
first  order  differential  equations.  Evaluating  the  limit  of  the  right 
hand  side  of  Eqn  (41)  leads  to 
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lira  .1406+2  -  420  v*  -  Sd*2^2 

xt-»0  (  + 

x 


lim 
x+  ■+  0 


(2806 


d6+  _  420|^  -  l06+v+2d5+ 


dx' 


dx"1 


ln»+^  +dv, 
-  105  v 

dx 


Combining  these  terms  with  the  left  hand  side  of  Eqn  (41)  yields  the 
following  initial  coefficients. 


5+  CH  v*  -  280) 


+2  4- 

C2  =  125^  vj  +  420 


sf  «5  -  K  > 


Similarly,  evaluating  the  limit  of  the  right  hand  side  of  Eqn  (42) 
leads  to 


lim 

x+-»0 


(240/Pr  -  Sd*2^ 


B 

x  ’  -*-0 


,  d(6+2v+) 

<_3^?  > 


lim 
"x  0 


(-65+v+ 


d5^ 

d? 
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+2  dv 
dx^ 


and  combining  these  therms  with  the  left  hand  side  of  Eqn  (42)  yields 
the  following  initial  coefficients  for  the  energy  equation. 


-  ►+  + 
74t’» 
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8  *+3  + 

-7  {i  vi 


Finally,  the  equation  of  motion,  Eqn  (41)  and  the  energy  equation, 
Eqn  (42)  can  be  solved  simultaneously  producing 


,(16+ 


c  /c  -  c  /c 

3i  2i  6i  5i 

C1  /C2  "  c4  /c5 

i  zi  *i  5i 


(47) 


W^i 


C,  /C,  -  C,  /C. 

3i  xi  6i  4i 

c  /c  -  c  /c 

i  i  i  i 


(48) 


With  the  information  supplied  from  Eqns  (41)  through  (48)  iterative, 
approximate  numerical  techniques  can  be  used  to  find  a  solution  to  the 
coupled  non-linear,  non-dimensionalized,  first  order  differential  equation 
of  motion  and  energy  equation,  Eqns  (41)  and  (42),  respectively.  A 
numerical  technique  is  outlined  in  Appendix  A. 
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V.  Results  and  Conclusions 

Laminar  Region 

Quantitative  data  of  the  laminar  boundary  layer  region  encompassing 
the  vertical  cylinder  are  extracted  from  Fig  11  using  section  III  and 
Havener's  computational  scheme  outlined  in  Ref  11.  The  wall  temperature 
is  52  C  and  the  ambient  temperature  is  17  C.  The  pressure  assumed 
constant  through  the  boundary  layer  is  approximately  0.1  mega-pascals. 
Using  the  film  temperature  to  evaluate  the  properties  of  air;  the 
Grashof-Prandtl  number  product,  GP,  is  approximately  3.24  X  10^  at 

8 

x  =>  .227  m.  Schlichting  (Ref  16)  theorized  laminar  flow  for  GP  <  10  . 

The  fringe  shift  data  extracted  from  Fig  11,  and  the  polynomial 
curve  fit  is  shown  in  Fig  12.  A  typical  temperature  profile  corresponding 
to  the  fringe  shift  data  at  x  =  .227  m  is  shown  in  Fig  13.  In  the  ambient, 
the  fringe,  which  is  a  line  of  constant  optical  path  length,  remains 
perpendicular  to  the  axis  of  the  cylinder  because  the  density  field 
is  constant  in  this  region.  However,  the  fringe  deviates  in  the 
boundary  layer  region  because  the  density  varies  through  the  boundary 
layer.  The  fringe  shift  provides  a  direct  calculation  of  this  deviation 
and  is  a  maximum  at  the  cylinder  wall  where  the  density  of  the  air  is  a 
minimum.  For  this  study,  the  fringe  shifts  are  calculated  along  a  fringe 
which  means  the  spatial  coordinates  of  these  shifts  do  not  share  a  common 
cross-sectional  plane.  This  approach  is  required  to  obtain  enough  data  to 
determine  a  good  fringe  shift  profile.  In  theory,  the  interferometric 
integral  equation  applies  to  planes  that  are  perpendicular  to  the  axis 
of  the  cylinder  so  the  shifts  should  be  planar.  The  non-alignment  of  the 
data  points  with  respect  to  the  x-coordinate  is  small  and  therefore  this 
approach  to  measuring  the  fringe  shifts  is  assumed  to  be  valid. 
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Pig.  11  Finite  Fringe  Interferogram  of  Laminar  Region 
Encompassing  Circular  Cylinder 


Fringe  Shift 


POLYNOMIAL  FIT 


Radial  Distance  y,  cm 


Fig.  12  Quantitative  Fringe  Shift  Data  at  Radial 
Distances  y 
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.382E-01 


Radial  Distance  y,  cm 


Pig.  13  Temperature  Profile  in  Radial  Direction  y  at 
x= . 227  n 


The  local  Nusselt  number  (NUx)cy^  and  convective  heat  transfer 
coefficient  h  at  x  =  .227  m  is  calculated  using  the  temperature  profile 
from  Fig  13  and  the  following  heat  balance  equation  where  k  =  2.888  X  10 
watt /cm- C 


q 

conduction 


-k  #) 
dy  w 


=  q  .  h(t  -  t  ) 

convection  w  » 


(49) 


2 

The  convective  heat  transfer  coefficient  is  6.73  watt/m  -C.  The  local 
Nusselt  number 


/dtv 


(Nu  >  ,  -  4*  -  (dy)wX 

ux  cyl  k  — z - 


v  - 


(50) 


is  approximately  56.9.  Using  the  wall  temperature  to  evaluate  the 
properties  of  air  since  (Nux)cyi  *s  proportional  to  the  temperature 
gradient  at  the  wall;  Gr^  =  3.648  X  10^  and  Pr  =  .7036.  A  theoretical 
value  for  (N^)^^  taken  from  the  curve  Pr  ■=  0.70  of  Fig  14  which 
is  reproduced  from  Sparrow's  and  Gregg's  numerical  results.  The  value 
of  £  at  x  =  .227  m  is  .4016  and  (N^)^-^/ (Nux)fp  =  1*18.  The  local 
Nusselt  number  for  the  flat  plate,  (Nu)fp  is  determined  using  Eckert's 
equations. 

(NUx)fp  =  0.508  Pr*4  (0.952  +  Pr)"*4  Gr^4 

(Nu^)fp  is  approximately  31.48  and  the  theoretical  (Nux)cy^  *s  37.15. 
The  relative  disagreement  is  34.7  percent  between  the  theoretical  and 
experimental  local  Nusselt  number. 

The  slope  of  the  numerically  derived  temperature  profile  will 


influence  the  experimental  Nusselt  number  and  the  magnitude  of 
disagreement  with  the  theoretical  Nusselt  number.  Error  induced  in 


1.5 


Fig.  14  Comparison  of  Local  Nusselt  Humber 
With  That  for  Flat  Plate 
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the  experimental  fringe  shift  profile  will  be  magnified  in  the  calculations 
for  a  local  Nusselt  number.  More  data  points  would  enhance  the  reliability 
of  the  experimental  results. 

Glassman  and  Orr  (Ref  8)  have  developed  a  method  to  automate  the 
reduction  of  data  from  an  interferogram.  Using  a  computer  and  an  auto¬ 
mated  optical  scanner,  the  fringes  of  an  interferogram  are  digitized  and 
local  values  of  the  fringe  shift  are  calculated  for  variable  mesh  speci¬ 
fications. 

Using  this  apparatus,  an  interferogram  of  the  laminar  region  (Fig  11) 
is  digitized  and  the  local  fringe  shift  values  calculated.  Subsequently, 
of  these  fringe  data,  a  contour  map  (Ref  1)  is  produced  to  obtain  lines 
of  constant  fringe  shift  throughout  the  field.  Following  this,  distribu¬ 
tions  of  fringe  shift  in  planes  normal  to  the  axis  of  the  cylinder  are 
determined.  From  here  on,  the  data  reduction  process  is  identical  to  that 
previously  referenced  in  this  section. 

The  technique  will  eliminate  error  caused  by  the  selection  of  non¬ 
planer  data  points  of  a  representative  cross-sectional  area.  The  spatial 
coordinates  of  the  relative  fringe  shifts  will  share  a  common  cross- 
sectional  plane.  Also,  the  number  of  data  points  is  increased  quite 
significantly. 


Transition  Region 

Of  the  three  regions  for  boundary  layer  analysis,  transition  is  the 
least  understood.  An  attempt  is  made  to  better  understand  the  transition 
region  both  qualitatively  and  quantitatively. 

Schlichting  uses  the  following  criteria  for  the  Grashof-Prandtl 
number  product  to  define  the  transition  region, 

108  <  GP  <  1010 

Figure  15  is  the  result  of  an  inter ferogr am  of  the  vertical  cylinder 

and  boundary  layer  for  a  wall  temperature  of  127  C  and  an  ambient 

temperature  of  19  C.  Using  the  film  temperature  to  evaluate  the  relevent 

properties  of  air  and  selecting  a  point  on  the  cylinder  corresponding  to 

the  middle  of  Fig  15,  the  Grashof-Prandtl  number  product  is  calculated 
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to  be  approximately  9.2  X  10  . 

It  is  interesting  to  note  the  eacistance  of  fringe  reversal  next  to 
the  cylinder  wall  in  Fig  15.  One  possible  explanation  for  the  fringe 
reversal  is  as  follows.  In  the  transition  region,  the  flow  begins  to 
become  unstable.  The  instability  at  the  onset  of  transition  is  caused 
by  small  perturbations.  The  perturbations  may  be  due  to  forces  or  momentum 
fluxes  contributed  by  the  ambient  surroundings.  These  perturbations  can 
strengthen  and  cause  the  boundary  layer  thickness  to  fluxuate.  At  the 
onset  of  transition  the  amplitude  of  the  fluxuatlon  is  small  and  increases 
in  magnitude  as  the  turbulent  region  Is  approached. 

For  laminar  flow,  the  gradients  for  temperature  and  density  are  the 
steepest  at  the  wall.  The  gradients  are  also  the  steepest  for  this  region 
as  compared  to  the  transition  and  turbulent  regions.  The  fluxuatlon  of 
the  boundary  layer  thickness  in  the  transition  region  influences  the 


Fig.  15  Finite  Fringe  Interferogram  of  Transition 
Region  Encompassing  Circular  Cylinder 
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steepness  of  the  gradients,  because  they  depend  on  the  boundary  layer 
thickness  for  convective  heat  exchange. 

The  fringes  are  lines  of  constant  optical  path  length.  The  optical 
path  length  is  proportional  to  the  index  of  refraction  integrated  over 
a  finite  path  length.  For  a  gas  the  index  of  refraction  is  proportional 
to  density.  So  the  optical  path  length  is  proportional  to  the  density 
which  is  inversely  proportional  to  temperature.  A  change  in  the  steepness 
of  the  density  and  temperature  gradients  influences  the  shape  of  the 
fringes  in  an  interferogram.  The  fringe  reversal  in  Fig  15  shows  this 
change  in  the  density  profile. 

The  fringe  reversal  can  also  be  observed  in  Fig  11  of  the  laminar 
boundary  layer  region.  These  qualitative  results  suggest  that  the 
fringe  reversal  may  be  evidence  of  the  phenomena  that  causes  laminar 
flow  to  transition  to  turbulent  flow. 

In  Fig  15,  the  fringe  reversal  appears  to  be  non-axisymmetric. 

This  may  be  caused  by  a  slight  tilt  in  the  cylinder  axis  with  respect 
to  true  vertical.  Havener  (Ref  11) ,  has  suggested  a  method  to  calculate 
the  error  of  a  non-axisymmetric  fringe  shift  compared  to  an  axisymmetric 
case.  Appendix  B  contains  information  on  this  error  analysis.  Using 
Havener's  method,  the  possible  error  incurred  here  is  3  percent. 

The  qualitative  data  obtained  in  this  investigation  supports  the 
theory  that  the  flow  is  periodic.  The  characteristics  of  the  transition 
flow  appear  to  represent  wave  phenomena. 

Figure  16  is  a  zeroth  fringe  interferogram  of  the  top  portion  of 
the  vertical  cylinder  and  free  convection  condition  as  shown  in  Fig  15. 

The  houndary  layer  is  non-laminar  as  indicated  by  the  irratic  inter¬ 
ference  pattern.  The  boundary  layer  thickness  varies  randomly  in  this 
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region  and  there  axe  locations  where  the  flow  appears  to  resume  or  Imitate 
a  laminar  flow,  figure  17  shows  the  same  conditions  as  fig  16  except  the 
Interference  pattern  is  now  a  finite  fringe  pattern.  This  random  Inter¬ 
ference  pattern  of  the  transition  boundary  layer  means  the  optical  path 
length  is  'varying  rapidly  in  this  region.  A  qualitative  assessment  of  this 
condition  is  that  the  density  is  also  varying  rapidly  but  not  randomly; 
while  the  local  pattern  is  irratic  and  indicates  the  existence  of  steep 
gradients,  there  also  appears  to  be  a  wave  structure  throughout  this 
section  of  the  transition  region.  Superposition  of  figs  16  and  17 
outlines  this  section  of  the  transition  boundary  layer  and  reveals 
qualitatively  this  wave  phenomena. 

for  this  investigation,  since  each  interferogxam  is  an  instantaneous 
record  of  the  flow,  transient  phenomena  is  difficult  to  identify  in  a 
single  exposure.  Accordingly,  data  of  the  wave  phenomena  can  be 
extracted  from  a  motion  picture  film  of  the  transition  boundary  layer. 
Schlieren  color  photographs  of  the  transition  region  reproduced  from 
motion  picture  film  exposed  at  500  frames  per  second  are  presented  in 
figs  18  through  34.  Every  tenth  frame  has  been  enlarged  and  printed 
for  qualitative  assessment  of  the  wave  structure;  these  coIot  enlargements 
definitely  suggest  that  this  transition  flow  is  periodic.  Each  dis¬ 
turbance  (a  wave)  that  propagates  along  the  cylinder  wall  appears  to 
he  followed  by  a  new  distrubance  that  follows  the  so  tailed  tracks 
qf  the  previous  disturbance.  This  periodic  phenomenon  is  visible  by 
comparing  the  first  and  eighth  photographs  where  the  new  disturbance 
can  be  seen  as  a  replacement  to  the  old  disturbance. 

As  a  disturbance  propagates*  along  the  cylinder  wall,  a  vortex  is 
Men  to  be  created  and  shed.  The  boundary  layer  flow  separates  into 

*  Observed  in  movie  of  boundary  layer  presented  with  thesis. 
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Schlieren  Color  Photographs  of  Transition  Region 
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Individual  mass  packets  and  some  o£  the  packets  form  into  the  vertex. 
The  vortices  are  periodic  and  start  and  terminate  at  a  prescribed  loca¬ 
tion.  It  is  the  author's  belief  that  the  vortices  are  caused  by  a 
momentum  flux  of  air  passing  from  the  ambient  to  replace  the  warmer 
air  that  has  dissipated  into  the  ambient  surroundings  due  to  the  shed 
vortices . 

The  special  fringe  reduction  technique  discussed  for  the  laminar 
region  could  be  extended  to  the  transition  region.  A  good  quantitative 
assessment  of  the  transition  region  could  be  achieved. 


Future  Recommendations 

In  a  method  developed  by  Havener  (Ref  2) ,  a  computer  is  used  to 
produce  finite  fringe  interferograms  from  theoretical  density  distributions 
of  a  field.  These  computational  interferograms  are  compared  to  actual 
test  interferogrmas  to  see  how  well  th-»  two  agree.  If  the  agreement 
is  good,  the  theoretical  density  field  may  be  assumed  to  be  the  same  as 
the  actual  field;  which  means  the  density  has  been  determined  without 
actually  reducing  experimental  interferometric  data.  For  the  case  studied 
in  this  research,  the  solution  to  the  set  of  first  order,  non-linear 
differential  equations  or  the  solutions  presented  by  Cebeci  and  Sparrow 
and  Gregg  could  be  used  to  produce  computational  interferograms. 

The  computational  interferograms  could  be  compared  to  the  test  inter¬ 
ferograms  to  confirm  the  numerical  solutions. 

Also  in  future  efforts,  the  contour' mapping  technique  could  be 
extended  to  analysing  the  turbulent  region.  If  a  density  function 
describing  the  turbulent  region  is  obtained.  Havener’s  method  could 
be  utilized  to  prepare  computational  interferograms.  The  information 
obtained  would  aid  in  the  understanding  of  turbulent  flow  phenomena. 
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Appendix  A 


for  Numerical  Solution  of  Governing 


System  of  Equations 


An  Iterative  technique  is  used  to  obtain  a  numerical  approximation 

4 

for  the  non-dimensional  boundary  layer  thickness.,  o  -  The  same  pro¬ 
cedure  is  also  used  to  obtain  an  approximation  for  the  non-dimensional 
velocity,  v+.  Both  of  these  parameters  are  functions  of  only  the 

Prandtl  number  Pr  and  the  non-dimensional  axial  coordinate,  x* . 

4.  4 

The  following  two  equations  are  used  to  obtain  6  and  v  at 
(x+  +  Ax+ ,  n) , 


*+  ,  d$  .  A  +  «.■+ 

6  -  ( — r  )  Ax  +  6 

dx  n-1  n-1 


dv  *  A  +  ,  + 

T+)  .  ^  +  Va 

dx  n-1 


(A— 1) 


(A-2) 


Here,  n-1  denotes  the  initial  values  at  x  =0. 

4  4 

After  values  for  £  and  v  are  known  at  the  new  station  n,  values 
n  n 

for  the  slopes  (d6*/dx* )  and  (fv+  /dx* )  are  cal  related-  Snowing  the 

n  ti 

values  of  6  ,  v  ,  d£  /dx  ,  and  dv  /dx  at  two  stations,  a  third  degree 
polynomial  numerical  approximation  is  used  to  calculate  better  values 
for  $*  and  v*  at  a  third  station  n+1  which  is  another  station  incremented 


from  n  by  Ax  .  The  general  form  of  the  polynomials  used  are 


a  +  bx_ 


+2  ,  j+3 

n+1  *  *»+! 


*’ + b'Ci -  ■=’=&  +  d'Ci 


<A-3) 


<A-4) 
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can  be 


Figure  A-l  displays  the  three  stations.  The  value  for 

evaluated  because  6+  6+  , ,  d6*/dx+,  and  d6+  n/dx+  .  are  known  at  two 

n,  n-1  n  n  n-1  n-1 

previous  stations.  Next,  the  coefficients  in  Eqn  (A-3)  are  evaluated 
from 


.  +  +  v  d<S  .  ,  d<5  .  0/»+  a-*-  % 

(x  -  x  ,)  — z  )  +  — z  )  -i  “2(6  -  o  ,) 

n  n-1  .  +  n  .  +  n-1  n  n-1' 

«-  dx  dx  J 


_ dx 

(xn  -  S-l> 


fdS+  ,  dS+.  |  ,  ,  +2  +2  . 


2<Xn  - 


Vl> 


b  -  -  2cx+  -  3dx+2 

.  +  n-1  n-1  n-1 

dx 


.+  u  +  +2  j 

a  m  6  .  -bx  ,  -cx  ,-dx  , 

n-1  n-1  n-1  n-1 

after  which  Eqn  (A-3)  is  used  to  calcualte  6+ at  station  n+1  or  6n+^. 

The  same  procedure  is  used  to  calculate  the  coefficients  in  Eqn  (A-4) 

■f  +  + 

and  v  at  station  n+l  is  also  calculated.  Next,  the  slopes  (d£  /dx  )n+^ 

+  + 

and  (dv  /dx  )  are  calculated  using  the  pair  of  coupled  equations, 
n+1 

Eqns  (43)  and  (44).  The  values  at  stations  n+1  and  n  are  subsequently 
used  to  calculate  new  values  at  the  next  station  which  is  a  Ax+  increment 
away.  Repeated  use  of  this  numerical  approach  is  carried  out  for  all 
&x+  increments  over  the  cylinder  length. 


The  coupled  equations ,  Eqns  (43)  and  (44)  are  a  set  of  non-linear 


first  order  differential  equations.  For  that  reason,  the  numerical 
approximations  for  S+,  v+,  d6+/dx+  and  dv+/dx+  are  highly  sensitive  to 
the  coefficients  in  Eqns  (A-3) ,  (A-4) ,  (43)  and  (44) .  The  coupled 
equations  are  also  singular  at  x+  =  0.  Further  understanding  of 
Eqns  (43)  and  (44)  is  required  to  obtain  correct  numerical  profiles  for 
6+  and  v+. 
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Appendix  B 


From  the  comparison  of  an  elliptic  density  field  of  eccentricity  e  to 
that  of  a  circular  axisymmetric  density  field.  Havener  (Ref  11),  developed 
the  following  equation  in  the  analysis  of  the  error  in  the  fringe  shift 
equation  (Eqn  1)  due  to  a  non- axisymmetric  density  field. 


S(y)e  is  the  fringe  shift  corresponding  to  the  ith  path  through  a  non- 
axisymmetric  elliptic  density  field  and  S(y)Q  is  the  fringe  shift  corres¬ 
ponding  to  the  ith  path  through  an  axisymmetric  circular  density  field. 
The  cross-section  of  the  axisymmetric  density  field  is  a  circle  that  is 
concentric  with  the  ellipse  of  the  cross-section  of  the  non-axisymmetric 
density  field.  The  circle  is  of  minimum  radius  required  to  encompass 
the  ellipse. 

If  an  elliptic  density  field  is  assumed  to  be  axisymmetric  the  error 


For  an  error  analysis  of  this  investigation  Fig  15  of  the  transition 
region  is  used.  Measuring  from  the  cylinder  axis  to  the  edge  of  the 
boundary  layer  (measured  at  the  reference  wire)  the  left  side  yields 
b  ■  1.75  units  and  for  the  right  side  a  ■  1.70  units.  Assuming  the 
cross-sectional  area  of  the  density  field  to  be  elliptic;  the  eccentricity 
of  the  ellipse  is  e  *  .24  and  the  error  induced  would  be  approximately 
3  percent. 
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A  Blight  tilt  In  the  cylinder  axis  will  not  necessarily  result  in  an 
elliptic,  non-axisymmetric  density  field.  The  cross-sectional  area  of 
the  density  field  fox  this  investigation  is  assumed  to  be  axisymmetric  so 
that  only  one  view  is  required  to  examine  the  density  field.  Multiple 
views  are  required  to  determine  the  shape  of  the  cross-sectional  area 
of  an  asymmetric  field.  Since  instantaneous  multiple  views  are  impossible 
to  obtain  for  this  study;  the  assumed  elliptic  shape  of  the  cross-sectional 
area  of  the  density  field  leads  to  an  approximate  3  percent  error  in  the 
data  reduction  of  the  fringe  shift.  In  the  scope  of  experimental  investi¬ 
gations,  this  error  is  tolerable. 
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Schlieren  methods  area  used  to  photograph  the  boundary  layer  to  provide 
a  qualitative  description  of  the  transition  region. 

An  approach  is  outlined  to  find  an  approximate  numerical  solution 
for  the  temperature  distribution  through  the  boundary  layer.  Integral 
equations  are  developed  for  a  differential  control  volume  within  the 
boundary  layer  after  which  assumed  velocity  and  temperature  distributions 
are  used  to  reduce  the  integral  equations  to  first  older  non-linear 
differential  equations.  The  differential  equations  are  non-dimensionali- 
zed  and  numeric al  techniques  are  vised  in  an  attempt  to  obtain  a  solution. 

Quantitative  data  from  the  laminar  region  resulted  in  a  discrepancy 
between  theoretical  and  experimental  evaluation  of  a  local  Ifusselt 
number.  Qualitative  observation  of  the  transition  region  reveals  a  fringe 
reversal  phenenona  that  oould  possibly  be  associated  with  the 
fundamental  mechanisms  of  transition. 
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